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, While multifractal spectra are convex, general field theoretic arguments 

Q\ [ show that power of field operators $ yield a concave spectrum of exponents 

as function of /. This is resolved by appropriate choice of operators to de- 
scribe multifractal moments. In a Lagrangian field theory of two mutually 
interacting species of fields <j),tf>, operators Off = ^ $ with traceless sym- 
metry give rise to multifractal spectra of harmonic diffusion near absorbing 
fractals when evaluated for zero component fields. 



> ■ I. INTRODUCTION 

l> . 

£N| , The concept of multifractality developed in the last decade has proven to 

be a powerful tool for analyzing systems with complex statistics which other- 
f — , ' wise appear to be intractable. It has found direct applications in a wide 

range of fields including turbulence, chaotic attractors, Laplacian growth phe- 
nomena etc. Here we generalize the idea of Cates and Witten by 
deriving the multifractal (MF) spectrum in the frames of a field theoretical 
^ \ (FT) formalism and make use of renormalization group (RG) methods. We 

relate the MF spectrum to the spectrum of scaling dimensions of a family of 
composite operators of Lagrangian 4 field theory. This gives an example of 
power of field operators whose scaling dimensions show the appropriate con- 



vexity for a MF spectrum, [11] while there is no need to include field gradients 
for this property. We calculate the MF spectrum to third order of perturba- 
tion theory using two complementary approaches: zero mass renormalization 
with successive e-expansion (see e.g. |12(| ) and massive renormalization group 
approach at fixed dimension, reproducing previous results obtained in 
lower order of perturbation theory for special cases. [B ,10 ] The resulting series 



are asymptotic. We take this into account and obtain numerical values only 
by careful resummation. 

We address a special case of a growth process controlled by a Laplacian 
field. The latter may describe a variety of phenomena depending on the 
interpretation of the field. For diffusion limited aggregation this field is given 
by the concentration of diffusing particles, in solidification processes it is given 
by the temperature field, in dielectric breakdown it is the electric potential, in 
viscous fingers formation it is the pressure. []5|Jl5[| In all mentioned phenomena 
the resulting structure appears to be of fractal nature and is characterized by 
appropriate fractal dimensions. [|J Its growth and spatial correlations lead to 
(non-trivial) spectra of multifractal dimensions. In general, the boundary 
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conditions determining the field will be given on the surface of the growing 
aggregate itself. It is this dynamic coupling that produces the rich structure of 
the phenomena and seems to make the general dynamical problem intractable. 

Here we study the simpler case when the fractal has been already formed 
and look for the distribution of the Laplacian field and its higher moments 
near the surface of the fractal. We will follow the diffusion picture, 

considering the aggregate as an absorbing fractal, "the absorber". The field 
p{f) gives the concentration of diffusing particles and vanishes on the surface 
of the absorber. More specifically we consider the Laplacian field p(f) in the 
vicinity of an absorbing (fractal) path, 'a polymer', or a junction of absorbing 
paths, 'a core of a polymer star'. In general we assume the ensemble of 
absorbers to be characterized by either random walk (RW) or self-avoiding 
walk (SAW) statistics. Multifractal scaling is found for the /-moments (pf (r)) 
of the field with respect to these ensembles. 

This formulation of the problem allows us to use the polymer picture 
and theory developed for polymer networks and stars [16-|l8[] and extended 
for copolymer stars. |19|,2Q| The theory is mapped to a Lagrangian (ft 4 field 



theory with several couplings [21-p3[] and higher order composite operators 
to describe star vertices. 

In section 2 we present the path integral solution of the Laplace equa- 
tion and relate it to polymer theory. Field theoretical representation and 
renormalization are discussed in section 01 together with the renormalization 
group flow and expressions for the exponents. In section IV we define the 
multifractal spectrum and give its series expansion in both renormalization 
group approaches adopted here. Section [V] is devoted to resummation of these 
asymptotic series and numerical results followed by some conclusions and an 



outlook in section VI 



II. PATH INTEGRAL SOLUTION OF THE LAPLACE EQUATION AND 

POLYMER ABSORBER MODEL 

In this section we describe the diffusion of particles in the vicinity of a 
polymer absorber by a "polymer" formalism. Let us formulate the problem 
first in terms of diffusion of particles in time. The probability of finding a 
randomly walking particle at point f\ at time t which started at point fb at 
time t = is given described by the following normalized path integral: 

G°(fo,n,t) = (5(^(0) -roW^it) -n))n (t). (1) 

Angle brackets in (jlj) denote the normalized integral 

_ /(...) eM-n (t))d{rW} 
^™') J exp(-n (t))d{rW} ' { ) 

which is performed with the Hamiltonian: 
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Integration in (|l|) is performed over all paths f^ 1 ^(r) with < r < t. Note that 
we have absorbed the diffusion constant into re-definition of time. G° (r*o , f\ , t) 
obeys the following differential equation: 

( A " + |-|) G ° (r »- r " () = ° 

with d the dimension of space. For finite d the random walker will visit 
any site after some finite time and we may assume a steady state limit for 
G°(ro,ri,t) for t = oo. In this case G°(fo,ri,t) will become independent of 
r*o and its limit and defines a field: 

P(ri) = lim / df Q G (r ,ri,t), 

here V is the system volume. The field p(r) obeys the Laplace equation: 

Ap(r) = 0. 

We introduce boundary conditions in such a way that the field p{r) equals 
to some constant p^, at r = oo and vanishes on the absorber. The absorber 
itself we describe by a path f^ 2 \s), < s < S^- 

Let us explain the solution of the Laplace equation (||) in the presence of 
an absorbing path r^ 2 \s). The boundary conditions are implemented by an 
avoidance interaction u\2 punishing any coincidence of the path rW of the 
RW and the path of the absorber. The correlation function of a random 
walk in the presence of an absorbing path r^ 1 (s) may then be written as 



G(ro,ri, Si) = (^(O) - r )5(r< 1 )(5 1 ) - n) 




where we have adopted the notation t = S\. 

We are interested in ensemble moments {p^ 2 (r*i)} of the field in the vicinity 
of the absorber, assuming an ensemble of RW or SAW absorbers. For the RW 
ensemble the average is performed with respect to the Hamiltonian Tio(S2), 
for the SAW ensemble an additional interaction has to be included. 

Here we choose a more general formulation, which allows us to describe 
the moments of the field in the vicinity of the core of an absorbing polymer 
star, or near the junction of f\ absorbing paths. 

The calculation of the /2 moment of (0) near the junction of /i absorb- 
ing paths will include the average over /2 random walks ending at r\ and 
the ensemble average over the configurations of the f\ absorbing paths with 
junction at T\. 

We are thus lead to consider the partition function of a star of walks which 
are in part mutually-avoiding. We will give this partition function here for 
the more general case of / walks of which fi random walks describe the field 
and /i walks correspond to the absorber. This situation describes the /2th 
moment of the flux of the field to the core of an absorbing star of /1 walks. 
We allow for additional avoidance interactions among these absorbing paths: 



1 TT -J') f 1 f f 

Z *f = ^o~([[ S (h (°) ~ ^o) exp | - - E ^afc / ds a / ds fe 

**/ i=l a,6=l J ° J ° 

5(r^(s a )-r< b Hs b )}) Kf{SaV 



here 



K f {s a } = J2n (s l 



i=i 



where ?io(Si) is given by @. stands for the partition function of star 
with zero interactions u a b = 0. The matrix u a b is given in the following form: 



Uab 



u\i if a, b < fi 

u° 12 ifa</i<6</i + /a 
or 6 < /i < a < /i + / 2 
else 



This corresponds to the partition function of co-polymer stars consisting of 
two species of chains |L£ . 2C ] with f\ chains of one species and f 2 chains of the 
other. is the interaction between absorbing paths. 



III. FIELD THEORY AND RENORMALIZATION 

As is well known, the polymer model may be mapped to the limit of m = 
of 0(m)-symmetrical Lagrangian field theory. |2(| To describe polymers and 
interacting random walks at the same time we adopt the formalism developed 



for multicomponent polymer solutions. [23] Its field theory is described by the 
following Lagrangian: 

/ 



a,a =1 

in general m-component theory 



C{<f>b,»b} = \ E fd d r (W* + (V^(r)) 1 

a=l 

1 1 r 

+ 4, E ««,«' J d d r^(r)^,(r). 



a=l 



fjL a is a chemical potential conjugated to the Gaussian surfaces S a in (|7|). 
Correlation functions in this theory are defined by averaging with the weight 
given by (0): 

((...)>Lc= / V[Mr)](-..)eM-£{<Pb^b}}\ m =o. 



here functional integration J T>[(j) a {r)} is defined in such a way that normal- 
ization is already included: (l)\c = 1 if all u aa > = 0. The limit m = in ( |l2|) 



can be understood as a certain rule to calculate the diagrams appearing in the 
perturbation theory expansions and can be easily checked diagrammatically. 

The partition function Z*f defined in (]7|) is mapped to the field theoretical 
correlation function Z*f via a Laplace transform in the Gaussian surfaces S a 
to conjugate chemical potentials ("mass variables") /j> a : 



POD 

Jo 



and 

/ 



Z*f{fJ-a} = ( / dr a J] 4>a(ro)<Pa(r a ))\c 



Our interest is in the scaling properties of these functions. Note that by 
( |j~4| ) these are governed by the spectrum of scaling dimensions of the composite 
operators Ila=i 4>a- To extract them we use renormalization group methods. 
[27, 28 1 Here we use the results of our previous approaches to the problem of 



co-polymer stars: massless renormalization group scheme with successive e- 
expansion (see e.g. |l2|]) and massive renormalization group approach at fixed 



dimension |13|] compiled in a pseudo-e expansion. j30[ On the basis of correla- 
tion functions it is standard to define vertex functions T„ corresponding to 
the couplings u a b as well as vertex functions rj^ a with insertion of composite 



operators Ila^a- Explicit expressions may be found in [19,gC]]. We define 



renormalization and introduce renormalized couplings g a b by: 

U a b = ^Z^ZfoZabgab. 

The renormalizing Z-factors are power series in g a b according to the following 
conditions: 

d 

Z0 a (g aa )—rW(u aa (g aa )) = 1 
Zab(gab)F^ab b (Uab(gab)) = f^9ab 

\i is a scale parameter (equal to the mass at which the massive scheme is 
evaluated and giving the scale of external momenta in the massless scheme). 

In order to renormalize the star vertex functions we introduce renormal- 
ization factors by: 

/ 



:E^ 2 )z^r^ a ( Uab (g ab )) = ^, 



^,1/2^ ry*S -p*/ 

a=l 

where Sn<f> a is the engineering dimension of the composite operator 



<W„ = /(f-l)+ 4- e 



The dependence of the renormalized couplings g a b and of renormalizing Z- 
factors on the scale parameter fi is expressed by the following relations: 



f^^dab = Pab(9a'b') 



d , w 



We are going to look on the situation of having two sets of walks of different 
species. In this case only three different couplings remain. We will refer to 
them as g\\, g 2 2, <7i2 = 521- The corresponding functions /?n, $22, $12 define a 
flow in the space of couplings. This renormalization group flow was discussed 
in [22,2*| . Its fixed points are determined by a set of equations: 



Pab(g*ab) = °, a, 6 = 1,2. 

In the space of the three couplings one finds [^] 8 fixed points corre- 
sponding to absence or presence of inter- and intra- species interaction. They 
are given in the table | where g* corresponds to the fixed point of the theory 
containing only 1 species, whereas gj. corresponds to the case of having only 
inter-species interactions, g^ describes a set of random walks interacting with 
another set of self-avoiding walks. The phenomenon we address in this article 
corresponds to the case of non-vanishing interaction between the two species 
of walks, while one set has no self-interaction. Thus we consider the two 
fixed points labeled G and U. The first corresponds to a set of random walks 
interacting with random walks of another species and thus describes absorp- 
tion on random walk absorbers, the second corresponds to a set of random 
walks interacting with another set of self-avoiding walks and thus describes 
absorption on SAW (polymer) absorbers. 

Having fx walks of the first species and /2 walks of second species we 
define the following exponents in the fixed points G, U: 

Vfif 2 = Vntaim = 922 = 0,512 = gh): 
r //i/ 2 = Vn^ign = g* ,922 = 0,512 = gh), 

which govern the scaling properties of the partition sum (|8|) . 

The scaling may be formulated in terms of the size R of the walks: We 
have to normalize the partition function by the number of configurations of 
the absorber given by -H*/^- For large R the resulting quantity scales like 

Z *hhl Z *W ~ R~ Xflh ■ for R = S u -»• oo 

Here \f 1 f 2 = f]f 1 f 2 ~ VfiO, u is the correlation length critical exponent of the 
walks: v = 1/2 for random walks and v ~ 0.588 for self-avoiding walks at 
d = 3. For the fixed point G we have r/^ = and v = 1/2 for all walks. 

The exponent A2, n corresponds to the nth moment of the flux onto the 
center segment of an absorbing linear chain. Considering the absorber to be 
either a random walk or a self-avoiding walk let us define the exponents: 

G _ J3 



^Rw{n) = A 2 = -r] 2 . 



^SAw{n) = A[L = V2n + V20 
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Previously [|lS,2C] we obtained the expressions for the exponents Vfafoi 
iH x * i n terms of e-expansion and pseudo-e expansion series in massless and 
massive renormalization group schemes. Whereas the first corresponds to 
collecting perturbation theory terms of the same powers of e = 4 — d, in the 
pseudo-e expansion series |3(| at each power of the pseudo-e parameter (r) 
one collects contributions from the dimension-dependent loop integrals of the 
same order. In the final results r = 1. Based on the expressions for the 
exponents r]f x f 2 , Vfcfo 1 1£ , 20 ] we find: 

n(n-l)e 2 n (n - 1) (-1 + n + 3 t(3)) e 3 
Aw(e) = en + ^ !! 



XsAw(z) 




69nC(3) 135n 2 C(3) 



+ ZTTTT^ E 



512 512 



\ I \ , / en 2 • £71 ■\2,\ 3loop 3 

arw{t) =ren+ [ — e n i t - — + e n% t ] r + X rw t°. 



3ren / enii 9 en 9 en i 1 eni 2 en\ 2 
W(T) = ^ + ^ + ^2 16 - + ^-T) T + 

\ 3loop 3 
A SAW T ■ 

Here £(3) — 1.202 is the Riemann zeta function, u are the loop inte- 
grals dependent on the space dimension d: at d = 3 i\ = 2/3, %i = —2/27. 
The expressions for the three-loop terms A^^ p , ^"aw m (l30| ) , (|3l|) are given 
elsewhere, [pll 



IV. MULTIFRACTAL SPECTRUM 

A widely used description for the MF spectrum is obtained from a Leg- 
endre transform, the spectral function f(a), of the analytically continued 
spectrum A(n) — nA(l) by 

/(a) = (A(n) — reA(l)) + na with a = — f-(X(n) — nX(l)) 

an 

For standard moments of a MF measure not including an ensemble average 
/(a) is called the spectral function. Here this notion is kept. The spec- 
tral function is widely used to characterize the multifractal nature of many 
processes. M In the standard approach the function /(a) defined for a multi- 
fractal measure on a set X gives for every a the fractal dimension of the subset 
of X for which the measure at scale I is characterized by t a with the Holder 
exponent a, in the limit £ — > 0. Due to this interpretation as a property of a 



(multifractal) measure, strict convexity conditions hold. The standard f(a) 
has the shape of a cap. While simple power of field operators $ will not 
generate such a spectrum, [11| the operators constructed here accord to this 



condition. 

Using the perturbation expansions for the A exponents given to third loop 



order both in e and r expansion in massless and massive renormalization (28) 
- ( |3l| ) and the relations for A(n) and the spectral function some algebra leads 
to the corresponding expansions for a n and f(a n ): 



n 3C(3) 3n 2 3nC(3) \ 3 

2 8~ ~8~ + ~4~ ' J ' 



orw(t) = 2+ (-| + eij + en-2enij)r 2 + a^ p r 3 



2 



f RW (r) = 2 + i^- - £ n 2 z, r 2 + /^V 



We) = 2- 1/4,+ (^"f> 2 

/ 149 69C(3) 21 n 81 n 2 135 <(3) 



+ - tttt - + — ' e 



/sAw(e) = 2 - l/4e + ( 



2048 512 512 512 256 
In 2 11 \ 



.3 



64 128/ 



]£ 2 



27 n 3 135n 2 C(3) 21 n 2 83 33C(3)\ 3 
256 + 512 1024 ~ 2048 + 256 r 



e ( e 9£nii 9n£ e%i eig 

" H 7- + ) T 



£ f £ HE Till 

a S Aw(r) = 2- -r +[---—— + 



16 4 16 



_2 



3/oop 3 
a SAW T 



Wr)-a- r+ (— + — + _ + -|-_i) T 2 + 

f Sloop 3 
•/ SAW 

Here £(3) — 1.202 is the Riemann zeta function, ij are the loop integrals 
dependent on the space dimension d: at d = 3 i% = 2/3, £2 = —2/27. 

Again ("(3) is the Riemann zeta function, i% and i% are the two-loop inte- 
grals and the explicit form of the three- loop contributions in (|35|), (|36[), (|39|), 



|40D is given elsewhere. [31] 
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V. RESUMMATION 



As is well known, the series of type (p8[)-(p|), (|33|) - (pCf), as they occur in 
field theory appear to be of asymptotic nature with zero radius of convergence. 
However, knowledge of the asymptotic behavior of the series as derived from 
the renormalization group theory allows us to evaluate these asymptotic series 



(see e.g. |29|). To this end several procedures are available depending on the 
additional information known for the series to be resummed. We assume 
our series to be of e-expansion character together with additional information 



for the case of the Lagrangian |23| ( |T0| ) we consider here. So we expect the 
following behavior of the kth. order perturbation theory term A\. for any of 
given above quantities: 



A k ~ k\k b (-a) 



k 



the constant a for the e-expansion of Lagrangian (f) field theory with one 
coupling was derived in |34| , |35|| : a = 3/8. For the unsymmetric fixed point, 
where two different couplings are present we use the value a = 27/64. [^3) We 
assume as well that the same properties hold also for the pseudo-e expansion in 
terms of r. With the above information in hand one can make use of the Borel 
summation technique improved by the conformal mapping procedure which 
up to now served a powerful tool in the field theory (see |2^] for example). 

Table [n| contains the results for the exponents \Rw{n) and XsAw(n) ob- 
tained in s and in pseudo-e expansion techniques from the corresponding 



values of exponents Vfi,f2 EH with the application of the resummation pro- 
cedure as described above. We have calculated the spectral function as a 
Legendre transform from the series of exponents received from the scaling of 
the moments of measure defined by diffusion. These moments were calculated 
as averages over all configurations of the absorber instead of performing a site 
average. Thus the interpretation of f(a) itself does not directly correspond 
to the picture developed above for the standard MF. All the same deriving 
f(a) in the same way from the spectrum of scaling exponents of moments of 
a measure the above discussed properties of a MF spectral function hold also 
for f(a) in this case. In particular max /(a) = tq gives the fractal dimension 
of the absorber and / is convex f"(a) < 0. Our numerical results for the spec- 
tral function are presented in Figs. |l]a,|l|b. They were obtained from the series 
for a n and f(a n ) as functions of n. We show the results of the resummation 
procedure described above applied to the series in both RG approaches. For 
comparison we also show the curve for direct summation of the e and r series 
to the 2nd order. In addition we have performed an analytical continuation 
of our series in form of [2/1] Pade approximants for the e 3 and r 3 series. It is 
obvious that direct summation of e 3 and t 3 series fails to converge and gives 
comparable values for a n , f{a n ) only for small values of n, i.e. near the max- 
imum of /(a) at n = 0. The symmetry of the Pade approximant holds only 
in the region shown and may be an artifact of the method. On the left wing, 
where it coincides with the resummed results the Pade approximant gives a 
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continuation which is compatible with the estimation for the minimal a value 
ctmin = d — 2. The Pade result is a m i n (e) = 1.333, a m i n (r) = 1.017 for the 
RW absorber and a m i n (e) = 1.250, a m i n (r) = 1.013 for the SAW absorber, 
which is calculated here only from 3rd order perturbation theory. 

Note that though the results obtained for a n and f(a n ) for a specific value 
of n differ in both approaches, the same curve f(a) is described with better 
coincidence for the left wing of the curves, corresponding to positive n. The 
resummation techniques we apply have proven to be a powerful tool already 
in the field theoretical approach to critical phenomena and have lead to high 
precision values for critical exponents. We hope that the application of these 
methods to calculations of MF phenomena allows to improve the reliability 
and comparability of the results. 

We currently work on the generalization of our present approach to find 
the series of spectral functions obtained for absorption at the core of a polymer 
with any number of arms. 

As can be extrapolated from the Pade approximant and as was shown also 
on the basis of high order approximations, || f(a) as it is defined here, will 
become negative near a m i n and a m ax- For this reason the identification of f(a) 
as the fractal dimension of some identifiable subset is not possible here. Also 
the extrapolation of the resummed data seems to indicate such a behavior. 
Note, however, that the perturbative approach, even in combination with 
resummation and analytical continuation is still not capable to give reliable 
results for high values of the expansion parameters. In particular this method 
is only good near the maximum of /(a). 

The possible negative values of the spectral function were discussed al- 
ready in and a physical interpretation of f(a) was given as a histogram of 
the measure /u plotted in logarithmic variables. In this interpretation nega- 
tive f(ct) indicate that the number of sites with a certain logarithmic measure 
a ~ In fj, decreases as the size of the absorber R increases. Thus for large R 
this number can only be defined by an ensemble average, as performed here. 

VI. CONCLUSIONS 

We have studied the characteristics of harmonic diffusion in the presence 
of a fractal absorber. We related the description of diffusing particles near 
absorbing paths to interacting walks. Following the model proposed by Cates 
and Witten || we used the polymer formalism to describe both the absorber 
and the random walks of diffusing particles. The ensemble of absorbers we 
considered to have random walk or polymer chain (self avoiding walk) statis- 
tics. The flux of diffusing particles onto such kind of the absorber and its 
higher moments generate a multifractal measure. |Q] The MF properties of 
this measure we described by the spectral function formalism. [|| 

We performed our calculations in the frames of a field theoretical ap- 
proach, relating our study to the study of scaling properties of composite field 
operators and defining their spectrum. Namely we show that the "copolymer 
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star" |0,|2^] - star operators Yl^ 1 lib V'o^fc i n a field theory with interactions 
and generate a spectrum of scaling dimensions which transforms to a 
MF spectrum with the appropriate convexity property. To calculate this spec- 
trum we have used the massless renormalization group scheme and massive 
renormalization group approach at fixed dimension. 

We give the explicit expressions for different kinds of exponents describing 
our problem (formulas (|§, @, @, (||) as well as for the 

spectral function (|34|), (|38|), (|36|), ( [40| ) in terms of power series in e = 4 — d 
and in pseudo-e expansion. All calculations were performed in the third order 
of perturbation theory. In particular, in the second order in e we recover 
previously obtained results. |J 

Special attention was payed to the fact that the series are asymptotic and 
have zero radius of convergence. We have used Pade approximants to obtain 
analytic continuation of the series under consideration for non-zero value if 
the expansion parameter. In addition we applied resummation techniques 
well approved in field theoretic calculations in order to obtain reliable infor- 
mation for the spectral function /(a), Holder exponent a and exponents A 
governing scaling behavior of averaged density moments of diffusing particles 
(see figs. [j]a, [TJb and table |ll|). While standard in field theoretical studies 
of critical phenomena, the resummation technique as to our knowledge was 
not applied in the theory of multifractals. We hope that our attempt will at- 
tract attention for this possibility in the context of other problems arising in 
the theory of multifractal measures as well as that the presumed accuracy of 
our results might evoke comparable efforts by numerical simulation. Further 
studies devoted to the set of spectra associated with diffusion near the core of 
absorbing stars with higher numbers of arms. We hope to gain more insight 
on the unsymmetric behavior of the spectral function and find the envelope 
of the family of spectra. 
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TABLE I. Fixed points for the interactions of a system of polymers of two species. 



G U Q U S G U U' S 

9ii g* g* g* g* 

922 g* g* g* g* 

912 g* G g*u gfr g* 
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TABLE II. Exponents \Rw{n) and XsAw(n) obtained in e and in pseudo-e expansion tech- 
niques. 



n ^rw{z) ^rw{t) ^saw{z) ^saw{t) 

1 0.99 0.99 0.71 0.71 

2 1.77 1.81 1.31 1.33 

3 2.45 2.53 1.86 1.92 

4 3.01 3.17 2.34 2.44 

5 3.51 3.75 2.78 2.94 

6 3.95 4.28 3.19 3.41 
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(Fig. la) 



(Fig. lb) 




FIG. 1. Spectral function f(a) for absorption on (a) RW and (b) SAW. Solid curves: 1 - [2/1] 
Pade approximant for e 3 results, 2- [2/1] Pade approximant for pseudo-e 3 results; dashed curves: 
3 - e 2 results without resummation, 4 - pseudo-e 2 results without resummation; stars - resummed 
e 3 results; boxes - resummed pseudo-e 3 results. 
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